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СТРУКТУРНЫЕ СВОЙСТВА 8САЬЕ-ГКЕЕ ГРАФА 
БАРАБАШИ- АЛЬБЕРТ 

Рассматривается зсаіе-ггее граф Барабаши-Альберт, описывающий большие се- 
тевые структуры типа Интернет. Вводится и исследуется фундаментальная мат- 
рица графа, точно описывающая свойства случайно выбираемых ребер. Выво- 
дится формула коэффициента кластеризации и предлагается метод его сепара- 
бельной настройки, не изменяющей распределения степени связности вершин. 
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1. Введение 

В последние годы резко возрос интерес к исследованию больших сетевых струк- 
тур (БСС), таких как Интернет, сложные социальные и биологические сети, инже- 
нерные коммуникации, транспортные сети мегаполисов, межмолекулярные взаимо- 
действия в различных средах и т. д. Необозримые размеры и сложность БСС обусло- 
вили рост разнообразия их стохастических моделей и становление теории случайных 
графов как самостоятельной научной дисциплины. 

Одним из наиболее широко известных видов случайных графов является граф 
Эрдеша-Реньи [1, 2]. Его генерируют на N вершинах: всякую их пару случайно, 
с вероятностью р, связывают ребром. Различные характеристики такого графа - 
коэффициент кластеризации, диаметр, вероятности появления тех или иных подгра- 
фов и др. выражены через его параметры N и р [3]. Локальная степень к связности 
его вершин имеет биномиальное распределение вероятностей, ее среднее значение 
(к) = р(ІѴ — 1). В дальнейшем нас будут интересовать бесконечные или очень боль- 
шие графы. При (к) = сопзі и N —>■ оо биноминальное распределение степени вершин 
графа Эрдеша-Реньи становится пуассоновским (отсюда второе название этого гра- 
фа — пуассоновский граф). Граф Эрдеша-Реньи является простейшим случайным 
графом: наличие связи между любыми двумя его вершинами определяется "чистой 
случайностью", независимо от наличия и конфигурации связей между другими вер- 
шинами. Для моделирования ряда реальных сетей корректно использовать его не 
удается. Это привело к появлению большого числа работ, посвященных исследова- 
нию других видов случайных графов [4-7]. 

Предложенная в 1999 г. в статье Альберта Барабаши и Реки Альберт [5] модель 
больших сетевых структур, формируемых по правилу предпочтительного связыва- 
ния (называемому еще правилом «богатый становится богаче»), завоевала за послед- 
ние десять лет огромную популярность среди исследователей больших сетей [8-11]. 
Эта модель представляет собой случайный динамический граф, выращиваемый из 
небольшого графа-затравки путем неограниченно повторяемых шагов добавления к 



нему новой вершины с т ребрами. Свободнвіе концві ребер каждой новой вершинві 
присоединяются преимущественно к вершинам, богатвім связями: вероятноств рі со- 
единения ребра с г-й вершиной графа пропорционалвна локалвной степени связности 
кі этой вершинві: 

(1) р { = кі/Е^ . 

С ростом графа Барабаши-Алвберт (графа БА) ряд его числоввіх характеристик 
сходится к стационарнвім значениям. Так, известна стационарная вероятноств 
того, что случайно ввібранная вершина графа БА имеет степенв связности к [12]: 

/ \ ^ 2т(т + 1) 

ѵ ; ^ к к(к + 1){к + 2)' 

Р. Алвберт и А. Барабаши назвали предложеннвій ими граф и его модификации 
зсаіе-ггее (безмасштабнвіми) графами, подчеркивая принадлежноств асимптотиче- 
ского при к — > оо распределения степеней связности вершин ос к~ 3 (2) (для моди- 
фикаций графа — распределения ос к~ а , а > 0) классу масштабно-инвариантнвіх 
распределений. В последовавших за [5] интенсивнвіх публикациях А. Барабаши и 
другие авторві убедителвно показвівают, что зсаіе-ггее графві адекватно описвівают 
многие сетеввіе структурві, изучаемвіе различнвіми прикладнвіми науками. С ростом 
числа подобнвіх исследований понятие зсаіе-ггее графов расширилосв и приблизилось 
к понятию графов, ввіращиваемых по правилу предпочтения. 

Актуальность развития теории зсаіе-ггее графов обусловливается широкой обла- 
стью их применения и значимостью получаемых на их основе результатов. К числу 
наиболее важных из них можно отнести оптимальные стратегии борьбы с инфекцией 
в информационных и социальных сетях, распространение методов теории критиче- 
ских явлений на задачи анализа устойчивости больших сетей и многие другие. 



2. Задачи, решаемые в статье 

В работах по использованию 8са1е-ггее графов для моделирования больших сетей 
используются, как правило, асимптотические распределения степени связности вер- 
шин. В ряде случаев, когда исследуется подмножество вершин, имеющих небольшую 
степень связности, необходимо учитывать точное, не асимптотическое распределение 
вероятностей. Для графа БА точное распределение (2) найдено в [12]. Для других 
графов, построенных по правилу предпочтения, точное распределение степеней вер- 
шин неизвестно. Для графа Б А и его модификаций неизвестны также точные распре- 
деления вероятностных параметров ребер. Поэтому многие структурные характери- 
стики этих графов, важные с прикладной точки зрения, остаются малоизученными. 

В статье устанавливаются точные распределения степени связности вершин для 
графов, выращиваемых по правилу предпочтения. Эти распределения включают 
формулу (2) как частный случай. Устанавливается точное распределение парамет- 
ров случайно выбранного ребра графа БА. Выводятся аналитические выражения для 
коэффициента кластеризации и других структурных характеристик графа. Предла- 
гаются методы калибровки графа по коэффициенту кластеризации. 



3. Распределение локальной степени связности вершин в графе БА 



Положим для простоты, что генерация зсаіе-ігее графа начинается с «затравки»— 
полного графа на 2т + 1 вершинах (т ^ 1). На каждом шаге і = 1,2,... генерации 
к имеющемуся ІѴ-вершинному графу добавляется одна вершина и связвівается т 
ребрами со случайно ввібраннвіми, в соответствии с вероятностями (1), вершинами 
имеющегося графа. Тогда на любом шаге і генерации графа минималвная степенв 
его вершин составляет к т і п = т, а средняя степенв (к) = 2т. 

Разобвём на шаге і множество всех вершин имеющегося графа на подмножества 
А к (слои) вершин, содержащие вершинві с одинаковой степенвю к (к — т, т+ 1, ...). 
Вероятноств ^ к слоя А к определим как вероятноств принадлежности этому слою 
случайно (равновероятно) ввібранной вершинві имеющегося графа: ц к = \А к \/М, где 
\А^\ — число вершин со степенвю к. Вероятноств Ц к того, что в бесконечном зсаіе- 
&ее графе случайно ввібранная вершина имеет степенв к, определяется как предел 
вероятности ^ к при і — > оо. 

При і — > оо вероятноств Р к связвівания нового ребра (одного из т ребер новой 
вершинві) с вершиной из слоя А к найдем из правила предпочтения (1): 

,„ч 0 ^іеА к к і \А к \-к А^ІѴ кЦ к 

і&А к 



Символом ~ будем обозначатв как сходимоств с нулевой относителвной погрешно- 
стью, так и сходимость по вероятности, в зависимости от контекста формулы. 

Для определения финальной вероятности Ц к выпишем уравнения баланса веро- 
ятностей при і — » оо. За один шаг генерации в граф добавляется одна вершина с т 
ребрами. В результате число \А т \ вершин слоя А т увеличивается на 1 (в этот слой 
попадает новая вершина) и одновременно уменьшается (т.к. из этого слоя уходят 
вершины, когда с ними связываются новые ребра) в среднем на тР т ~ т 2 ■ Я т /(к) 
вершин. Приравнивая выражения для доли вершин слоя А т в числе всех вершин 
графа до и после шага генерации, получаем уравнение баланса: 

\А т \ \А т \ + 1-т 2 д т /{к) 



N N + 1 

из которого, умножая его на ІѴ(ІѴ + 1), приводя подобные члены и применяя для 
избавления от переменной ІѴ асимптотическое равенство |А т |/ІѴ = <5 т , находим: 

(4) Я т 



{к) + т? 



Аналогично выписывается и решается уравнение баланса для С} к при к > т. 
Число \А к \ вершин слоя А к за шаг генерации в среднем (из-за прихода вершин из слоя 
А к _і) возрастает на тР к _і ~ тік — 1)Як-і/ (к) и уменьшается на тР к ~ тк(^) к / (к) . 
Соответствующее уравнение баланса 

\А к \ _ \А к \ + т(к - 1)Як-і/ (к) - гпкОк[ (к) 
N ~ N+1 



дает для финальных вероятностей соотношение: 

^ ~ тп ■ (к — 1) 

(5) к>т + 1 . 

Из рекуррентного определения (4), (5) последовательности нетрудно с учетом 
равенства (к) = 2т получить (посредством индукции) явное выражение (2). 
Примечания. 

1) Применяемая техника записи и решения уравнений баланса является сжатой 
формой записи для <2%(ІѴ) уравнения в конечных разностях о^(ІѴ) = ^(А+І) — 5&(ІѴ) 
и перехода к его стационарному решению Як = д%(оо), определяемому условием 
ІѴ -> сх) , гід А (ІѴ) = 0. 

2) В [5] и ряде других работ для /с — > оо используется оценка ~ 2т 2 к~ 3 в 
форме плотности /(5) = 2т 2 з~ 3 , з ^ т. Эта смещенная оценка вероятностей дает 
несмещенную оценку средней степени вершины графа: (к) = ] з>т з/(з)оІз = 2т. 

3) В программном генераторе графа БА выбор вершины новым ребром, соот- 
ветствующий вероятности (1), целесообразно выполнять в два этапа: вначале по 
распределению (3) следует разыграть номер к слоя, а затем выбрать любую вер- 
шину этого слоя равновероятно. Эксперименты с генерацией графов размером от 
500 до 50000 вершин показали, что такой генератор работает на порядок быстрее 
генератора системы моделирования АпуЬо§іс (кйр: //у/у/у/, х^гек. сот/) и в десятки 
раз быстрее генератора распространенной библиотеки для работы с графами ЛШС 
(Ьіір:/ /^ип§.8оигсеюг§е.пеі/). Этот метод ускорения можно применять в генераторах 
других видов графов, формируемых по правилу предпочтения. 

4. Обобщенная формула распределения локальной степени связности 

Если при построении графа задавать вероятность выбора г-й вершины новым реб- 
ром в виде Рі = /(/сі) / ^2 і /(^)) применяя функцию предпочтения / ^ 0 (любую), то, 
как вытекает из уравнений баланса, аналогичных приведенным в п. 3, распределение 
вероятностей принимает вид: 

а т ■ /(/с — 1) 

(6) Ят=— 77 — V' Як = Чк-і— тттѵ 

а + т ■ / (т) а + т ■ / (к) 

где (если такое распределение существует) параметр а = (/{к)) можно найти 
численными методами из уравнений ^2 к>т Як = 1 или ^2 к>т ^Як = 2т, или 
Т1ік>т ${к)Як = а. Например, т.к. логарифмическая / = Іод(к) и инверсная ло- 
гарифмическая / = 1/1од(к) функции предпочтения индуцируют ряд (6), отвеча- 
ющий уравнению ^2 к>т Як — 1 независимо от основания логарифма и значений 
а > 0, т > 0 (вещественных), то для определения параметра а = (/{к)) уравнение 
2~2к>т Як = 1 непригодно, и нужно использовать одно из двух оставшихся. 

Разные / позволяют получать разные (не только степенные) распределения С^ к . 

5. Задача о свойствах случайно выбранного ребра графа Б А 



Вообще говоря, асимптотически-степенное распределение степени связности вер- 
шин, положенное в основу наименования 8са1е-ггее графов, слабо определяет их 



структурные свойства. Это можно видеть из приведенных на рис. 1 примеров гра- 
фов, которые все имеют одно и то же степенное распределение степени связности, 
но определенно не все являются графами БА. Слева изображен граф БА, по центру 
- зсаіе-ггее дерево, и справа — зсаіе-ггее «творог», состоящий из бесконечного числа 
компонент связности. Степень связности вершин каждого из трех графов распреде- 
лена по закону (2) при т = 2. 



Рис. 1. Разнообразие структур при фиксированном распределении степени связности 

Написанные нами генераторы зсаіе-ггее дерева и зсаіе-ггее «творога» реализуют 
методы, отличные от правила предпочтительного связывания. Из приведенных при- 
меров видно, что структурные свойства графа БА определяются методом его генера- 
ции, и что распределение степени связности является только необходимым, но не до- 
статочным признаком адекватности описания конкретной исследуемой сети графом 
Б А. Чтобы углубить статистическое согласование генерируемых графов с моделиру- 
емыми сетями, используют различные модификации правила предпочтения. Наряду 
с распределением степени связности [12] при согласовании учитываются такие струк- 
турные характеристики, как коэффициент кластеризации графа [13], диаметр графа 
или число подграфов определённой структуры [14]. Подобные структурные характе- 
ристики оцениваются приближенными асимптотическими методами или путем тру- 
доемкого статистического моделирования, поскольку для точного их расчета необхо- 
димо учитывать совместное распределение (неизвестное) степеней вершин на концах 
случайно выбранного ребра. Ниже задача о свойствах случайно выбранного ребра 
графа БА решается путем вывода рекуррентной формулы, определяющей матрицу 
О, совместного распределения концевых степеней связности дуг в ориентированной 
версии графа БА. Матрица О, позволяет поставить нахождение конфигурационных 
характеристик графа БА на систематическую основу, и в этом смысле является его 
фундаментальной матрицей. 



Введем ориентированную версию графа БА, чтобы зафиксировать различие 
свойств у двух концов любого его ребра. Это различие возникает в момент добавле- 
ния ребра в граф вместе с новой вершиной в силу того, что только один конец ребра 
выбирает вершину, имеющую высокую степень связности. В ориентированной версии 
граф будем выращивать добавлениями новой вершины с т дугами, начала которых 
инцидентны этой новой вершине, а концы присоединяются к имеющимся вершинам 
графа, выбираемым по правилу предпочтения (1). Такая замена ребер дугами не из- 
меняет графа БА, но несколько упрощает его описание в ходе его анализа. В этой 




6. Ориентированная версия графа БА 



используемой далее нотации граф БА представляет собой случайный динамический 
ориентированный граф С = С(і). 

Слой Ак графа С(і) определяется по-прежнему как множество вершин, имеющих 
степень связности к. Вершины слоя А т (граничного) по построению инцидентны 
только исходящим из них дугам. Случайно выбранная дуга (с. в. д.) характеризуется 
степенью (1,к), где I — степень вершины, из которой дуга исходит (начальная степень 
дуги), к — степень вершины, в которую дуга заходит (концевая степень дуги). Рас- 
пределение вероятностей степени с.в.д. опишем матрицей О. = вероятностей 
того, что дуга имеет начальную степень I и концевую степень к (1,к ^ т). 

Определим туннель В(1,к) как множество дуг, начало у которых инцидентно 
вершинам слоя Аі, а конец — вершинам слоя А^. Туннели В(т,к) с параметром 
/ = т будем называть входными, с параметром I ^ т + 1 — внутренними. 

Всякий слой Ак на очередном шаге эволюции графа пополняется новой верши- 
ной тогда, когда новая дуга выбирает вершину предшествующего слоя А^-і (если 
он не пустой), и, наоборот, теряет вершину, когда новая дуга выбирает вершину это- 
го слоя А^, т.к. выбранная вершина переходит в слой А^+х- Новая вершина всегда 
добавляется в граничный слой А т . 

Число дуг, исходящих из вершины любого слоя А}., равно т, число заходящих 
в нее дуг равно к — т. Число дуг туннеля В(1, к), заходящих в случайно (равнове- 
роятно) выбранную вершину слоя А^, лежит между 0 и к — т, и в среднем равно 
\В(1, к)\/\Ак\, где \В(1, к)\ — число дуг туннеля В(1,к). 

При і — » оо число вершин N — > оо, |А&|//Ѵ ~ С^^, \В(1, к)\/(тМ) ~ Яі,к, г Д е 
тіѴ — число всех дуг в графе, (^і^ — вероятность того, что с.в.д. принадлежит 
туннелю В(1, к). С учетом этих свойств графа С(і) выпишем уравнения баланса для 
определения (^і^ (1,к ^ т). Техника составления и решения уравнений баланса для 
вероятностей дуг аналогична используемой в п. 3. 



7. Вероятности входных туннелей 

7.1. При добавлении новой дуги (скажем, первой из т дуг новой вершины) ее 
конец попадает в слой А^ (к ^ т + 2), согласно (3), (2), с вероятностью Р& ~ = 

= ]]г^щ^2), где (к) = 2т — средняя степень вершины графа. Выбранная дугой 
вершина вследствие повышения степени переходит в слой Аь +1 и выводит из входного 

г>( 1 \ \В(т,к)\ \В(т,к)\ 

туннеля В[т, к) те дуги, которые в нее заходили — в среднем 1 ^ и ~ 1 фщ^ ДУ Г - 
С учетом вероятности Р^ выбора слоя А& среднее число дуг, выводимых новой дугой 

тэ< 1 \ (га+1) \В(т,к)\ 

из туннеля В(т, к), сходится к (А .| 1)(А .| 2) • щ^ 1 - 

7.2. Число дуг в туннеле В(т, к) уменьшается и тогда, когда конец новой дуги 

попадает в слой А т , т.к. при этом вершина, в которую вошла дуга, переходит в слой 

А т+1 и выводит из входного туннеля В(т, к) те его дуги, которые из нее исходили. 

Поскольку новая дуга выбирает вершину в слое А т случайно, то среднее число ИСХО- 
тз/ і \ \В(т,к)\ т\В(т,к)\/(тЮ тС} т к 

дящих из этой вершины дуг туннеля В[т, к) равно 1 гд | = ^ А гщ - ~ ^ ' . С 

учетом вероятности Р т ~ ( т ^^т+2) этого события и соотношения (2) среднее число 
дуг, таким образом выводимых из туннеля В(т, к), сходится к Р т ■ т ® тЛ ~ т Ят,к ^ 

7.3. Найдем теперь среднее число дуг, добавляемых в туннель В(т, к) новой дугой 
графа. С вероятностью Р^-і ~ ^~7&р^ = ыи+і) коне Д новой дуги попадает в неко- 
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торую вершину слоя А^-і, к ^ т + 2. Эта вершина, переходя в слой А^, переводит в 
состав туннеля В(т, к) в среднем т^~у заходивших в нее дуг туннеля В(тп, к — 1) 
и саму добавленную новую дугу. С учетом вероятности Рк-і этого события новая 
дуга добавляет в туннель В(тп, к) в среднем |^±^ + ^л^ГГ^) Д У Г - 

7.4. Приравнивая формально выражения долей дуг туннеля В(т, к) в числе всех 
дуг графа С(і) до и после добавления новой дуги, получаем уравнение: 



(7) 



\В{т,к)\ _ \ В ( т > к )\ + ц0)( 1 + 



(т+1) ( 1 _,_ \В(т,к-Х)\ \ _ (т+1) \В{т,к)\ _ тОт^ 
МЯ к -г ) (к+\){к+2) ЫС! к 2 



гаіѴ тіѴ + 1 

Умножая его на тИ^тИ + 1), приводя подобные члены, и вновь деля на тіѴ, имеем: 
|Б(т,/с)| (т+1) ^ , |Б(т,А;-1)|^ (т + 1) |Б(т,/с)| т^ т , 4 



(тЫ) к(к + 1)\ ЛГд к _! У (&+1)(А; + 2) ІѴд 

тЦ т ^ (т + 1) / т|5(т, /с — 1)| ^ (т + 1) т|і?(т,А;)| 



Дг(ЛЧ-І) V 7 (Л + 1)(А; + 2) тіѴф* ' 

(т + 2) (т + 1) Л тфт.Е-Л (т + 1) т(^ ті/г 

^ т '" 2 = А;(А; + іД + ^ (Л+1)(А; + 2) ' 

/ а ч п (т + 1) / , гпОт^к-Л ( т + 2 т(т+1) ^ -1 

Ш/с>\ гл 2т(т+1) /"ѵ 2т(т+1) 

дставляя в (8) Ц к = щ^щ^щ и = ^зщійт) , получаем рекуррентную 

формулу, определяющую всю последовательность €} т ^ через <З т , т +і: 

2(т+ 1) (А; — 1) „ 

Начальный элемент <5 т , т +і последовательности найдем из специального урав- 
нения баланса для входного туннеля В(т,т + 1). Для него, т.е. при к — т + 1, 
повторяются рассуждения, приведенные выше для туннелей В(т, к) при /с ^ т + 2, 
с единственной поправкой. Она состоит в том, что при выборе концом новой дуги 
слоя Ак_і (в данном случае — слоя А т ) в туннель В(т, к) = В(т,т+1) добавляется 
лишь одна эта новая дуга. Соответствующий член уравнения (7) и выражения (8), 
учитывающий добавку \В(т, к — І^ДіѴС^.!), при к = т + 1 исчезает, и получаемое 
уравнение баланса решается так: 



(ю) я 



т,т+1 



(т + 2)(2т + 3)' 

Формулы (10), (9) определяют всю последовательность (5т,*;> к ^ т + 1. 



8. Вероятности внутренних туннелей 

Внутренние туннели В(1,к), 1,к ^ т + 1 начинаются и заканчиваются во внутрен- 
них слоях графа. Уравнения баланса для них составляются на основе рассуждений, 
подобных рассуждениям о входных туннелях, с учётом следующих двух отличий. 



Первое отличие состоит в том, что при добавлении к графу новой дуги наряду с 
рассмотренными прежде возможностями пополнения туннеля В(1,к) появляется еще 
одна возможность. А именно, когда конец новой дуги попадает в вершину слоя Аі-і, 
эта вершина переходит в слой А\ и переносит в него начала дуг, концы которых уже 
лежат в слое А^. Эти дуги тоже пополняют туннель В(1, к). Среднее число таких дуг 
определяется выражением ^^ 1 '^ . С учетом вероятности Рі-\ новая дуга добавляет 

таким способом в туннель В(1, к) в среднем Р;_і ^^ 1 '^ ~ ^ж^Т^ ДУ Г- Второе 

отличие сводится к тому, что сама новая дуга никогда не добавляется во внутрен- 
ний туннель В(1,к). Уравнения баланса для внутреннего туннеля получаем, внося 
в (7) изменения, учитывающие два отмеченных отличия, и заменяя фиксированный 
индекс т индексом I: 

(п) Щ^ = , 1 ,-Ым+І т , + Т щ !: к ~ 1) ^ 



тЫ) (тіѴ) + 1 ѵ ' У| к{к+1) ІѴд^ 

т + 1)\В(1 -1,к)\ (ш + І) \В(1,к)\ (т + 1) \В(1,к)\ 



1(1 + 1) ІѴд^ (к + 1)(к + 2) (/ + 1)(/ + 2) ЫС^г 

Решая (И) с помощью приемов, использованных при решении (8), находим: 

т(т + 1) Оі^к-і т(т + 1) Оі-і^ т(т + 1) т(т + 1) (^і^ 

Чі ' к ~ к{к + і) д*_х /(/ + і) д,_! [к + і){к + 2)я к (і + і)(і + 2) я г ' 

Наконец, подставляя сюда для С^к,Як-і,ЯьЯі-і, их выражения в виде (2), получаем: 

(12) Яі,к = І + к + 2 ' ' = т 1 ' т + 2 '--- • 

Вместе с решениями (10), (9) и с учетом того, что в граничный слой А т дуги не 
заходят, решение (12) рекуррентно определяет вероятности всех туннелей орграфа 
С = С(оо) и тем самым матрицу С} =|| Яі,к || следующим образом: 



(із) Яі, к = < 



'0, I ^ т, к — т, 

(т+2)(2т+3) ' I = ГП, к = Ш + 1, 

2(т+1) . (к-1) п ] — т ѣ > т -1- 9 

к(к+1)(к+т+2) ^ (к+т+2)Ч1,к-1, ' — "Ч ^ "Г 

(й-і)<Эг^-і+(г-і)<Эг-і,й 



(г+/с+2) 



I ^ т + 1, к ^ т + 1. 



Расчет матрицы выполняется, начиная с заполнения крайнего левого столбца 
(столбца к = т) нулями и расчета элементов верхней строки (строки / = т): сначала 
по формуле (13) вычисляется ее второй слева элемент Ц т ^ т +\, затем, через него - 
следующие вправо элементы Ц т ,к-> к ^ т + 2. После этого построчно, слева направо, 
рассчитываются элементы следующих строк матрицы. 

На рис. 2 показан начальный фрагмент матрицы при т = 2. В табл. 1 приве- 
дены оценки вероятностей (^і^, полученные при т = 2 с помощью имитационного 
моделирования (ИМ) графа С(оо) (число N вершин в сгенерированном графе равно 
100 тыс.). 
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СП 

•и 

© 



1 
Й 

1 

со 

1 

-I I 

со 

1 
1 

V- 
3 1 

с 

1 

>< 

1 

Й 

^ 

зс 5 

сі 
сз & 

Гп зс 



0.071429 
0.017857 
0.005952 
0.002381 
0.001082 
0.000541 



0.064236 
0.020238 
0.007857 
0.003506 
0.001732 
0.000924 



0.050794 
0.018254 
0.007835 
0.003731 
0.001987 
0.001116 



0.039683 
0.015512 
0.007143 
0.003652 
0.002014 
0.001177 



0.031385 
0.01 2987 
0.006294 
0.003363 
0.001927 
0.001164 



0.025253 
0.01 0878 
0.005478 
0.003030 
0.001790 
0.001111 



0.020668 
0.009169 
0.004755 
0.002704 
0.001638 
0.001040 



10_ 

0.017183 
0.007792 
0.004136 
0.002405 
0.001 487 
0.000962 



Рис. 2. Начальный фрагмент матрицы О., рассчитанной по формуле (13) при т = 2 



Таблица 1 



Статистические оценки вероятностей (^і ъ при т = 2 





2 


3 


4 


5 


6 


7 


8 


2 


0 


0.071 


0.065 


0.050 


0.040 


0.032 


0.025 


3 


0 


0.018 


0.021 


0.018 


0.016 


0.013 


0.011 


4 


0 


0.007 


0.009 


0.008 


0.008 


0.007 


0.005 



Соответствие рассчитанных значений (^і^ (при т = 2 и других значениях т) их 
оценкам, получаемым путем ИМ, подтверждает правильность формулы (13). 

9. Маргинальные распределения начальной и концевой степеней дуги 

Исходящий из слоя А\ пучок туннелей В(1, *) определим как множество всех дуг, 
исходящих из вершин слоя А\. Множество всех дуг, заходящих в вершины слоя А^, 
образует заходящий пучок В(*,к). Вероятности С^і* и принадлежности с. в. д. 
пучку В(1,*) и, соответственно, пучку В(*,к) определяют маргинальные распреде- 
ления начальной и концевой степеней дуги соответственно, и могут вычисляться по 

формулам = Т,к>тЯі,к , Я*,к = Еі>тФ,*- 

Более простое выражение для (З^* найдем, замечая, что число дуг в пучке В(1, *) 
— дуг, исходящих из слоя Аі, — в среднем равно <27ІѴт (т.к. из всякой вершины графа 
исходит т дуг, а в слое А\ находится в среднем эдТѴ вершин). Отсюда, с учетом (2), 

\В(1, *)| д ь Ит 2то(га + 1) 



(14) 



Нт 



пт — — 

лг-»оо І\т 



Яг 



I ^ т . 



Щ + Щ + 2У 
В среднем начальная степень с. в. д. равна 2т. 

Случайный выбор дуги из множества дуг графа С(оо) эквивалентен случайному 
выбору вершины из множества его вершин, и последующему случайному выбору 
исходящей из нее дуги. Попутно заметим, что, поскольку средняя степень вершины 
(к) = 2т, то в среднем вершина имеет т исходящих и т заходящих дуг. 

Для нахождения ф*^ учтем, что число дуг в туннеле В(*, к) — дуг, заходящих в 
слой Ак, — в среднем равно (к — т)^кN, т.к. во всякую вершину этого слоя заходит 
к — т дуг, а в слое А& находится в среднем д&ІѴ вершин. Отсюда имеем: 

(к — т)д^ (к — т) 2{к — т){т + 1) 



(і5) д м 



Нт 

ІѴ— >оо 



Як 



к > т . 



тИ т к(к + 1)(к + 2) 

Из (15) следует, что в среднем концевая степень с. в. д. бесконечна. 
Начальная и концевая степени с.в.д. зависимы, т.к. произведение их маргиналь- 
ных распределений (14) и (15) не тождественно их совместному распределению (13). 
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10. Характеристики случайно выбранного ребра графа БА 

Если дуги орграфа С заменить ребрами, то он превратится в традиционный 
неориентированный граф БА. Свойства случайно выбранного ребра (с. в. р.) графа 
БА можно определить через свойства с. в. д. орграфа С, принимая во внимание, что 
фиксированный (любой) конец такого ребра с равной вероятностью является нача- 
лом или концом дуги. Так, вероятность того, что концевые степени ребра равны 
I и к (для случайно выбранного и для противоположного его конца соответственно) , 
можно определить формулой: 

(іб) ѳ^ = ѳ м = (д г ,, + д м )/2. 

Маргинальные вероятности концевых степеней с. в. р. в соответствии с (14), (15) со- 

Ѳгл I гл і гл \ /о (к—т)(т+1)+т(т+1) (т+1) 
к,* = ®*,к = [Чк,* + Ч*,к)/2 = кІк+і)(к+2) = (к+і)(к+2) и совпадают с 

финальной вероятностью Рк выбора слоя А% ребром новой вершины. С учётом (16) 

и (13) отсюда следует, что концевые степени с. в. р. статистически зависимы. 



11. Изменение среднего числа слоев при построении графа 

Для простоты числом слоев будем называть текущий наибольший номер слоя 
К = к тах . Слой А к с номером К = к тах может после появления в нем первой 
вершины быть выбран с вероятностью К \Ак\/^^к^ = А/ (2т А) одной из т дуг 
новой добавляемой в граф вершины. На каждый шаг выращивания графа малая 
(при больших А) вероятность такого выбора возрастает в т раз (по числу дуг) и 
составляет Рк = К/(2Ы) ос К/И. Т.к. этот выбор означает увеличение числа слоев 
К на единицу, то в среднем ЛК = Рк, т.е. сіК ос К/Ы, откуда следует, что К ос А^ 1 / 2 . 

При ИМ графа всегда возникает ситуация, когда выбранная новой дугой вершина, 
являясь единственной вершиной последнего слоя А%, уходит в следующий новый 
слой Ак+і и оставляет предыдущий слой Ак пустым. С ходом времени эта ситуация 
закрепляется: пустыми становятся несколько предпоследних слоев, а последний слой 
содержит ровно одну вершину. При этом вероятность выбора новой дугой вершины 
последнего слоя становится выше вероятности пополнения предпоследнего слоя. 

Оценки для А'ср(А), вычисленные путем ИМ при т = 2, представлены на рис. 3 
маркерами. Непрерывной кривой изображен график функции К ср = 3.05 А 1 / 2 . 



12. Коэффициент кластеризации графа Б А 

12.1. Коэффициент кластеризации С определяется как утроенное отношение сред- 
него числа пд треугольников в графе к среднему числу пу вилок (путей длины 2): 
С = Зпд/пу [15, 9]. Для приложений графа БА актуальна задача нахождения ана- 
литической формулы этого важного показателя [14]. Применение асимптотической 
пропорции К ос А 1 / 2 позволяет найти решение этой задачи в следующем виде: 

/ \ ^ (т-1)(1пА) 2 

(17) т>2. 

На рис. 4 изображены графики линеаризованных зависимостей от А коэффици- 
ентов кластеризации С, рассчитанных с помощью ИМ графа БА (здесь у = ЫС, 
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х = 1п[(1п 7Ѵ) 2 /ІѴ]). Форма графиков (параллельные прямые с единичным наклоном) 
подтверждает установленную теоретически функциональную связь (17). 

12.2. При т = 1 треугольники в графе не появляются, т.е. С(ІѴ) = 0 . Для вывода 
асимптотической формулы (17) коэффициента кластеризации рассмотрим вначале 
случай т — 2. 

Пара дуг новой вершины, добавляемой в граф, может увеличить число треуголь- 
ников в нем на 1 в том случае, когда дуги этой пары выбирают две смежные вершины 
графа. Если одна новая дуга выбирает слой Аі, а другая — слой , то вероятность 
р^ к получения при этом треугольника равна отношению числа дуг, соединяющих 
слои А[ и Аь, к числу всех возможных выборов пары вершин: 

д Г {|в(*,А:)| + |в(м)|}/{|Л| х \М}~™Щ<э 1)к + д кг1 )/(<э& к №), к ф Х ^ 

Рі > к ~{ В(1, /) |/ { ІМхМЬИ | ^ 2тМЯ ц /\Я 1 Н{Я 1 Н - 1)], к = І. 

В любом случае отсюда с ростом N имеем: 

(18) Р&~™(Яі,к + Яы)/(ЯіЯкЮ> 1,к>т. 

С учетом (18) и вероятности того, что концы двух новых дуг попадают в слои 
А\ и А к , новая вершина добавляет в граф в среднем ~ Рі,к т {Яі,к + Як,і) /{ЯіЯи^) 
таких треугольников, две вершины которых выбраны в слоях А\ и А\~. Общее число 
треугольников, добавляемых одной новой вершиной, в среднем составляет: 

1=т ^—^к=1 

где К (И) ~ сіѴ 1 ^ 



(20) 



1,к 



2Р ь Р к1 I ф к, 
РіРкі I — к, 



/о\ г> (т+1) г) (т+1) 

и, согласно (3), Р ~ ( Д 1)( ^ 2) > П ~ д^щщ - 
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С учетом (20) можно изменить диапазон значений индекса к в формуле (19) и 
переписать ее следующим образом: 



1=т ' ' к=т 



т 



(т + 1) 



(т + 1) 



X 



і=т 2-^к= т [(і + і)(і + 2)(к + 1)(к + 2) 
Щ + 1)(/ + 2) к(к + 1)(к + 2) (Ои + д м ) 
2т(т + 1) 2т(т + 1) 



х 



ІѴ 



\_ у-ВД у-ВД іНОі^ + Як,і) _ 1 у-ВД ^Х(ТѴ) 



,к- 



Далее, используя асимптотики ~ т(т+1)1 2 к 2 (см. (П. 6), прил.1), и К ~ сіѴ 1 / 2 
(см. п. 11), находим: 



(Іп А Ш) ~ > > 



к=т 



ІкѲ 



1 



І.к 



2тИ 



К{Ы) Г К(И) 



Ѳі ^ік ■ &Ык 



— / / т(т + 1)Г 2 к- 2 1к-Шк= У ' / А; -1 / Г 1 ■ Шк 



2тИ 
(т+1 



2ІѴ 
(т + 1) 



т 



т 



(т + 1) 
2ІѴ 



(т + 1)/, сіѴ 1 / 2 
ш 



2ІѴ 



' с ІпіѴ 

1п 1 

т 2 



где с' = 1п — (слагаемое С ± ^ в асимптотическом представлении приращения 

<іпд сохраняем в целях получения хорошего приближения (23) для конечных ЛГ). 
Суммируя приращения с/пд по времени і эволюции графа, находим среднее число 
треугольников пд в графе с N вершинами при т = 2: 



{т + 1) 
2ІѴ 



„/ 1п ІѴ 



(ІпіѴ) 2 с , ЛТ 
- — + 1п — 1п N 
4 тп 



2ІѴ \ т 

т+1)(1піѴ) 2 ДпіѴ 
8ІѴ +С 1\Г' 



(21) п А = ^гіпдф ~ / 

4=1 ^ 1 



(т + 1)(1пі) 2 



(т + 1)(1пА0 3 ,(1п ІѴѴ 

= ~ ^ 1" с — ^ — 



12.3. При т > 2 малая при больших N вероятность йпд формирования тре- 
угольника новой вершиной графа возрастает по сравнению со случаем т = 2 в 
С"т — тп(т — 1)/2 раз (по числу пар добавляемых дуг). При этом вероятностью об- 
разования сразу двух и более треугольников можно пренебречь в силу меньшего ее 
порядка малости. Соответственно, в раз возрастает среднее число треугольников 
в графе. Отсюда, учитывая (21), получаем: 



(22) 



т(т — 1) 



Пд 



(т + 1)(1піѴ) 3 ДІпЛО 



+ с 



и полученное выражение охватывает рассмотренные выше случаи т = 1 и т = 2 

12.4. Асимптотическое выражение для среднего числа пу вилок можно получить 
исходя из того, что каждая вершина в слое Ль является центром С\ 



к(к-1) 

К 2 ; вилок, 
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и, т.о., 



Ѵ- Сѵ ^ л^2 Г^^лг^и лг 2т(ш + 1) - 1) 



'*=™ * У т " к * У™ *(* + і)(* + 2) 

7Ѵт(т+1) / - 1 — ; — -^А; ~ ІѴт(т+1) / -<йк ~ Ытіт+І) ІпГсѴіѴ) = 

Ут (А + ІДА + 2) ] т к 

/77 ( 777 _|_ 1 ) 

= ЛГш(т + 1) [1п(с) + 0.5 1п ІѴ)] 1 ^ 1п ІѴ . 

12.5. Отсюда с учетом (22) находим коэффициент кластеризации графа БА: 



Зт(т— 1) 



(т+1) (ІпГѴ) 3 /{ІпТѴ^ 



т-1)(1піѴ) 2 ЬІѴ 



(23) С = — 1 , -г ^ = - — 1 + (т - 1)с т - л . . 

ѵ ; п у т(г ^ +1) ІѴ1піѴ 8 ѵ ІѴ 

Формула (17) непосредственно вытекает из (23). Коэффициент с т для конкретных т 
можно вычислять путем ИМ. Так, при т = 2, 3, . . . , 8 его значения приблизительно 
равны 0.5, 0.35, 0.275, 0.18, 0.12, 0.10 и 0.07 соответственно. С ростом т коэффициент 
с т убывает по экспоненте. 

13. Сепарабельная калибровка графа по коэффициенту кластеризации 

На практике при ИМ больших сетевых структур выясняется, что согласованный 
с ними по распределению степеней связности вершин граф БА нередко существенно 
отличается от них по коэффициенту кластеризации. Например, у сети Интернет этот 
коэффициент в несколько раз выше, чем у соответствующего графа Б А. Для согла- 
сования графа БА с моделируемой сетью по коэффициенту кластеризации можно 
модифицировать алгоритм генерации графа способом, который легко усматривается 
в ходе вывода формулы (23) и при т = 2 состоит в следующем. 

Вначале каждое из двух новых ребер выбирает себе слой по распределению (3) в 
соответствии с ускоренным алгоритмом генерации графа (см. примечание 3 в п. 3). 
Затем с вероятностью 1 — Л для выбранных слоев А\ и выполняется процедура 
случайного (равновероятного) выбора пары і, ] не смежных вершин г 6 А\ и ^ € 
(если такие имеются), с которыми и соединяются новые ребра (если такой пары не 
смежных вершин нет, то выбирается равновероятно любая пара г € Аі, ] € А^). Но с 
вероятностью Л выбирается случайно пара смежных вершин г Е Аі, ^ Е и новые 
ребра связываются с ними, образуя треугольник (если пары смежных вершин нет, 
то выбирается случайно любая пара г € Аі, ] € А^ и треугольник не формируется). 

Наименьшее значение 0 регулируемого коэффициента кластеризации С Т достига- 
ется при Л = 0. При Л = 1 значение С г достигает максимума С г = С тах , который 
при т = 2 можно найти следующим образом. При больших N и А = 1 каждая новая 
вершина добавляет треугольник во всех случаях, кроме одного — когда обе ее дуги 

выбирают слой А]~ = А т = А 2 . Т.к. вероятность этого случая Р| = ^ д^щ^щ ) = 

( ( ш +Т)(т+2) ) = іб' ТО с Р е Д нее число треугольников при Л = 1 составляет пд ~ Щ^. 
И, поскольку среднего числа вилок пу ~ т ( г ^+ 1 ) дг } п дг = ЗЛПпіѴ метод калибров- 
ки не изменяет, то С Т = С тах ~ ^ІмЬѵ = ІбТІЬѵ = пЙѵ^' г Д е С ~ коэффициент 
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кластеризации не калибруемого графа. Т.о., при 10 3 ^ ІѴ ^ 10 6 предлагаемый метод 
позволяет устанавливать любое значение коэффициента кластеризации в диапазоне 
от 0 до (20 3000) х С. 

Возможны различные способы обобщения метода на случай т > 2. В общем слу- 
чае параметр Л (0 ^ А ^ 1) может быть задан в виде функции Л = А(ІѴ). Метод 
можно применять и в генераторах других видов графов, формируемых по правилу 
предпочтения. Из построения метода видно, что управление коэффициентом класте- 
ризации в ходе генерации графа не влияет на распределение степеней вершин. 

14. Заключение 

Современная теория случайных графов характеризуется интенсивной разработ- 
кой и обоснованием математических моделей и методов, которые применяются для 
исследования больших сетевых структур, изучаемых прикладными науками. Одной 
из наиболее известных и удачных моделей является стохастический зсаіе-ггее граф 
Барабаши- Альберт (граф БА), описывающий структуры типа Интернет. Получен- 
ные в статье новые результаты углубляют и уточняют существующие представления 
о структурных свойствах графа БА. К числу основных теоретических результатов 
относится точное совместное распределение вероятностей для параметров случайно 
выбираемых ребер. Для его нахождения вводится ориентированная версия графа 
Б А, обеспечивающая возможность эффективно использовать технику составления и 
решения уравнений баланса для вероятностных параметров случайно выбираемых 
дуг. На основе такого подхода получена рекуррентная формула, точно задающая 
совместное распределение вероятностей начальной и концевой степеней дуги, а это, 
в свою очередь, позволило определить точное совместное распределение вероятно- 
стей для концевых степеней ребер неориентированного графа БА и асимптотические 
свойства этого распределения. 

К частным результатам, имеющим самостоятельное значение для дальнейшего 
исследования и применения зсаіе-ггее графов и графа БА, относятся: 

- метод генерации графов по обобщенному правилу предпочтения, использую- 
щему произвольную функцию предпочтения / ^ 0, и найденное в общем виде точное 
распределение локальной степени связности вершин в графах, генерируемых этим 
методом; 

- ускоренный метод генерации графов предпочтения, позволяющий сократить 
время их генерации по сравнению с существующими методами на порядок и более; 

- маргинальные распределения концевых степеней дуг/ребер графа Б А; 

- асимптотически степенной с показателем 0.5 закон роста максимальной степени 
вершин в процессе эволюции графа; 

- аналитические выражения коэффициента кластеризации и метод его сепара- 
бельной настройки, не влияющей на распределение локальных степеней вершин и на 
совместное распределение концевых степеней ребер графа. 
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ПРИЛОЖЕНИЕ 1 
Асимптотические характеристики дуг и ребер 



1. Общая формула элементов первой строки (строки I = т) матрицві О., опреде- 
ляемой рекурсивным соотношением (13), имеет вид: 

(П1) О 2 2(ш + 2)(ш + 3)...(2ш + 2) 

1 ^ Чт ' к ~ к{к + 1) к(к + 1)...(к + т + 2) ' ^ 

Действительно при к = т решение (П.1) дает 

2 2(ттг + 2)(т + 3)...(2т + 2) 



Яг 



т(т + 1) т(т + 1)...(т + т + 2) 



при к = т +1 из (П.1) получаем фт.т+і = ( т+2 )(2т+з) ' Т - е - П Р И этих ^ Р ешение (П.1) 
уравнениям (13) удовлетворяет. Остается проверить, что при любом к ^ т + 2 для 
(П.1) выполняется равенство: ф т , й = щ^^щ + щ^щЯт,к-і, т.е. что 

2 2(т + 2)(т + 3)...(2т + 2) 



*;(*; + 1) &(А; + 1)...(А; + т + 2) 

2(т + 1) (к-1) 



к(к + 1)(к + т + 2) (к + т + 2) 



2(т + 2)(т + 3)...(2т + 2) 



(Л; — 1)Л: (А; - !)&...(& + т + 1) 



Но это равенство является тождеством. Таким образом, решение (П.1) доказано. 
2. Общая формула элементов второго столбца (к = т + 1) матрицы О. имеет вид: 

2 ш ( ш + і)( ш + 2 )...(2ш + 2) 

(771 + 2) і(/ + 1)...(/ + 771 + 3) 

Действительно, при / = т из (П. 2) имеем (5 т)ГП +і = ( т+2 )(2т+з) ' ЧТО соответств У ет 
(13). При всяком I ^ т + 1 , согласно (13), решение (П. 2) должно удовлетворять 
условию Яі,к = ( ^ 1)<Э ''^^ 1)(3 '~ 1Л , где к = ш + 1, <2г,ь-і = <5г,ш = 0, т.е. должно 

/\ (/ — ™_і_і 

выполняться равенство Яі >т +і = - — (/ +т+ ^ , или, иначе, равенство: 

2 т(т + 1)...(2т + 2) (/-1) 2 т(т + 1)...(2т + 2) 



(т + 2) /(/ + 1)...(/ + т + 3) (/ + т + 3) (т + 2) (I - 1)1. ..(I + т + 2) ' 

Но данное равенство очевидно. Решение (П. 2) доказано. 

3. Общая формула элементов строки I матрицы О. имеет вид: 

2 777 ( 777 -Ь 1 ) 

где = 0(/с- т " 3 ) > 0 и при к = т значение Щ(к) = Д,(т) = 2/[/(/ + 1)]. 

Действительно, при / = т (первая строка матрицы О.) найденное выше точное 
решение (П.1) имеет вид (П.З). Покажем, что если для некоторой строки I — 1 ^ т 
общая формула имеет вид (П.З), то и для строки I она также имеет вид (П.З). Для 
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этого заметим, что при к = т (первый элемент строки I), формула (П.З) безусловно 
представляет верное решение С}і <т = 0. Остается убедитвся, что если представле- 
ние (П.З) имеет место для (^)і^-і, то оно справедливо и для (^і^- Подставляя в (13) 
Я^-і = щ"і"к-1 )к - Чк - 1), <2г_і Л . = { ^щк+і) ~ имеем: 

1 



(І + к + 2) 



2т(т + 1)(к-1) 2т(т + !)(/-!) 

I. щ + і)(к-і)к ~ ік ~ 1) - Яі{к (і-і)ік(к + і) ~ (/ ~ 1} • Кі ~ і{к \ 

1 \2т(т + 1)(1 + к + 2) 



(І + к + 2) [ 1(1 + 1)к(к + 1) 
2 п(т + 1) (*-1) д, « - 1) ^ = 2 т (т + 1) _^ 



/(/ + 1)Л(А; + 1) (І + к + 2) 14 у (/ + & + 2) 1 1Ѵ у /(/ + 1)А;(А; + 1 

где 



и, т.о., для строки / также справедливо представление (П.З). Свойство Яі(к) 

) - 1- 



0(Аг т - 3 ) ввітекает из (П.4) с учётом того, что Пг-і(к) ~ 0(&- т - 3 ), ~ 1 - ІЛА 



иД,(А-1)~ВД. 

4. Элементві столбца /с матрицві имеют следующее асимптотическое представ- 
ление (ввівод которого для краткости опускаем; оно легко проверяется численнвіми 
расчетами) : 

(П.5) Я ик = 8 к {1) ~ С%_ г С{т)1- т - \ к>т+1, 

где С к п _ 1 — биномиалвнвій коэффициент, С(т) = 2т(т+1) - [(т + 3)(т + 4)...(2га + 2)]. 
При к = т + 1 формула (П.5) непосредственно согласуется с точной формулой (П. 2). 

5. Асимптотика свойств с. в. р. при увеличении его концеввіх степеней I и к опре- 
деляется формулами (П.4), (П.5) и (16) в виде: 

777 I 777 -1- 1 1 

(П.6) <ді, к = Ѳ к ,і = (Яі, к + Я к ,і)/2 



1(1 + 1)к(к + іу 

Ѳі >к ~ т(т + 1)Г 2 к~ 2 . 

Сопоставляя (П.6) с маргиналвнвіми вероятностями Ѳ^* и Ѳ*^ (см. п. 10), заключа- 
ем, что с ростом т совместное асимптотическое по I и к распределение сходится к 
распределению независимвіх случайнвіх величин. 
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